Diagrammatic Analysis of the Unitary Group for Double Barrier 
Ballistic Cavities: Equivalence with Circuit Theory 



A. L. R. Barbosa and A. M. S. Macedo 

Departamento de Fisica, Laboratorio de Fisica Teorica e Computacional, 
Universidade Federal de Pernambuco 50670-901 Recife, PE, Brazil 

We derive a set of coupled non-linear algebraic equations for the asymptotics of the Poisson 
kernel distribution describing the statistical properties of a two-terminal double-barrier chaotic 
billiard (or ballistic quantum dot). The equations are calculated from a diagrammatic technique for 
performing averages over the unitary group, proposed by Brouwer and Beenakker [J. Math. Phys. 
37, 4904 (1996)]. We give strong analytical evidences that these equations are equivalent to a much 
simpler polynomial equation calculated from a recent extension of Nazarov's circuit theory [A. M. 
S. Macedo, Phys. Rev. B 66, 033306 (2002)]. These results offer interesting perspectives for further 
developments in the field via the direct conversion of one approach into the other. 

PACS numbers: 05.45. Mt, 73.23.-b, 73.63.Kv 



I. INTRODUCTION 

Random-matrix theory has proved to be a powerful 
tool for describing generic features of quantum chaotic 
systems^*?. For closed systems it offers an accurate sta- 
tistical description of both energy levels and wavefunc- 
tions characteristics with overwhelming numerical and 
experimental confirmation. In this case, the central hy- 
pothesis is to replace the Hamiltonian of the system by a 
matrix with random independent gaussian entries obey- 
ing certain fixed exact symmetries, such as time-reversal 
and spin-rotation. The resulting ensembles are known 
in the literature as the Gaussian Ensembles. For open 
systems, such as ballistic chaotic cavities and disordered 
wires attached to external leads, the major role is shifted 
from the Hamiltonian to the appropriate matrix associ- 
ated with the description of electron transport, viz. the 
scattering matrix or the transfer matrix. Here, the jus- 
tification of the specific ensemble is a subtle procedure, 
because, unlike the previous case, one needs to take into 
account subtle correlations between the entries of the ma- 
trix induced by flux conservation and quantum diffusion 
(for disordered wires). 

For ballistic chaotic cavities, a general procedure to 
obtain such an ensemble was developed by Mello and co- 
workers^. It consists of evoking the maximum informa- 
tion entropy principle along with certain generic assump- 
tions, such as analyticity, unitarity and some specific 
symmetries that are exactly preserved in the presence of 
chaotic dynamics. The resulting S-matrix distribution, 
P{S), turns out to be a multidimensional generalization 
of the Poisson kernel, usually found in two-dimensional 
electrostatics. Its general form reads 



P(S)(x \det(l-&S)\ 



(1) 



where (3 € {1,2,4} is a parameter identifying distinct 
symmetry classes^, M = 2N is the total number of open 
scattering channels (N for each waveguide) and S is a 



sub-unitary matrix, i.e. the eigenvalues of SS' may be 
less than unity. For a two-terminal system, consisting of 
a ballistic chaotic cavity or quantum dot coupled to two 
waveguides, the scattering matrix can be conveniently 
written as& 



5' 



r t 



(2) 



where t, t' and r, r' are respectively transmission and re- 
flection matrices. The presence of barriers of arbitrary 
transparencies at the interface between the chaotic cav- 
ity and the waveguides can be accounted for by specifying 
the average (or optical part) of the scattering matrix^. 



n o 

r 2 



(3) 



where r\ e r 2 are reflection matrices for barriers 1 and 2 
respectively. 

In mesoscopic physics, one is usually concerned with 
the description of a well defined measurement, such as 
the full counting statistics (FCS) of a two-terminal sys- 
tem, whose cumulant generating function is given by the 
Levitov-Lesovik formula 5 . 



$(A) = -M Trln(l + (e iA - l)tt f ), 



(4) 



where Mq = eVTo /h ^> 1 is the number of attempts to 
transmit an electron during the observation time To, t is 
the transmission matrix and V is the voltage. Physical 
observables can be obtained from the series expansion 



fc=i 



k\ 



where 



d(i\) 



A=0 



(5) 



(6) 
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are the irreducible cumulants of the FCS. Some well 
known examples are the dimensionless conductance and 
shot-noise power, which are given respectively by 

g = gi /M = Tr(«t), (7) 

which is the Landauer formula, and 

p = ? 2 /M =Tr[ttt(l-ftt)]. (8) 

The third cumulant has also attracted some recent in- 
terest, including an experimental observation in tunnel 
junctions 6 . It is given by 

k = q 3 /M = Tr[ttt(l - - 2tt f )]. (9) 

The above expressions for the cumulants of the FCS 
are sample specific, which means for a chaotic cavity that 
an average over the Poisson kernel distribution is neces- 
sary in order to make comparisons with real experimental 
data. Such calculations can in principle be performed ex- 
actly for arbitrary values of M and (3 using the method 
of supersymmetryiS and some very general results are in 
fact already available in the literatur o^ 1 ? . Albeit pow- 
erful and general, this method can become very cum- 
bersome and sometimes unwieldy. Fortunately, if one is 
interested only in the scmiclassical regime, where one ne- 
glects quantum interference contributions, defined math- 
ematically by the asymptotic condition M 3> 1, much 
simpler alternative techniques exist. 

Two well-known methods to deal with the semiclas- 
sical regime are the diagrammatic analysis of the uni- 
tary group, presented by Brouwer and Beenakker— and 
Nazarov's circuit theory 12 . In its original form circuit 
theory does not allow for a direct comparison with the di- 
agrammatic approach in the entire range of parameters, 
i.e. for arbitrary values of the barrier's transparencies, 
because of the intrinsic difficulty to determine the aver- 
age pseudo-current- voltage characteristics of an arbitrary 
circuit element. This is related to the well known prob- 
lem of breakdown of semiclassical transport equations 
close to boundaries and interfaces, where pure quantum 
effects become dominant. This difficulty was recently 
removed by a novel systematic treatment presented in 
Ref. ITjj . in which circuit theory is combined with the su- 
persymmetric non-linear sigma model yielding a power- 
ful technique with a very convenient algebraic structure. 
This extended version of circuit theory raises the natural 
question as to whether its information content coincides 
with that of the diagrammatic method in all ranges of 
the input parameters. This is a highly non-trivial ques- 
tion, because the semiclassical concatenation principle, 
which is the basis of circuit theory, does not have a di- 
rect representation in the diagrammatic formulation. In 
particular, there is no obvious way how to extend this 
concatenation principle to account for quantum correc- 
tions, such as weak- localization. We remark that the 
semiclassical concatenation principle is directly related 
to the independence of the leading asymptotic contribu- 
tion on the symmetry parameter [3. By contrast, the 



weak localization correction is strongly (3 dependent 
and, in particular, vanishes for f3 = 2, i.e. for systems 
with broken time-reversal symmetry. Interestingly, for 
the particular case of symmetric barriers a direct com- 
parison between both approaches was presented in Ref. 
Hi and complete agreement was found. 

Motivated by this result, we present in this work a de- 
tailed comparison between the diagrammatic approach 
and the extended version of circuit theory for the case of 
asymmetric barriers with arbitrary transparencies, thus 
completing the analysis of Ref. [l3| and providing strong 
evidence for the full equivalence between these semiclas- 
sical techniques. This result is particularly relevant in 
practical applications because of the great algebraic ad- 
vantage of the calculations in circuit theory (in our case 
study the problem is reduced to a polynomial equation of 
fourth order) in comparison with alternative approaches. 
In fact, we believe that circuit theory equations reach 
the ultimate irreducible form in terms of simplifying the 
description. Furthermore, when combined with the scal- 
ing theory presented in Ref. for the balistic-diffusive 
crossover in phase-coherent metallic conductors, in which 
the present results enters as an initial condition, we end 
up with a very powerful formalism for performing con- 
crete calculations in realistic conductors. Applications 
might include spintronics 14 and normal-superconducting 
hybrid systems 1 ^. 

II. THE DIAGRAMMATIC TECHNIQUE 

In this section we present a detailed account of the 
application of the diagrammatic method to the calcula- 
tion of the semiclassical limit of a double barrier chaotic 
ballistic cavity. It contains the central result of this pa- 
per and complements previous analysis by Brouwer and 
Beenakker— . 

We start by noting that all cumulants of the FCS can 
be generically written in the form 

A = Tr[a{ttf)}, (10) 

where a(x) is an arbitrary smooth function. It means 
that the ensemble averages of transport observables can 
be fully described by the following density function 

p(r) =<TW(r-ttt)), (11) 

so that 

(A) = f dTa{r)p{r). (12) 
Jo 

The problem is thus reduced to the calculation of p(r) in 
the asymptotic regime M = 2N ^S> 1 . 

In the diagrammatic formalism, one starts by intro- 
ducing the generating functions 

F 1 (z) = (C 1 (z-S^C 2 SC 1 )- 1 ), (13) 
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and 



where 



F 2 (z) = (Ciiz-SdS^d)- 1 ), 



(14) 



where S is the M X M scattering matrix defined in Eq. 
(J2J and the auxiliary matrices Gi and C 2 are defined as 



Ci 




1 



JY 



C 2 



In 




(15) 



in which 1jv is the N x N unit matrix. It can be easily 
verified that 

P (t) = — Im[/i(T + iO+)] = --Im[/ 2 (r + iO+)], (16) 

7T 7T 

where fi(z) = Tr(Fi(z)),i = 1,2. 

The next step is to define the matrix generating func- 
tion 



F{z) -(f 2 (z) 



(17) 



which can be written as a sum of two terms that can 
be averaged separately, F(z) = (2z)~ 1 (F + (z) + F-(zj). 
Each term is given by 



F a {z) = C + aA(l - ^ a Gl)- x ^ a B 



(18) 



where A = CL, B = aTC, with a = z~ 1/2 and the 
following matrices have been defined 



C = 



d 

c 2 



L = 



L 
it I ' 



and 



T \ fl-( R Q 

1 - I o T t J > W - I fit 



(19) 



(20) 



The submatrices L, T and R describe transmission and 
reflection coefficients of the barriers and are related to S 
of Eq. via the condition that 



U = 



S L 
T R 



(21) 



be unitary. The remaining matrices in Eq. (|18|) are the 

"free propagator" 



G°=R + aafCL 



(22) 



and the "self-energy" matrix, E CT . They satisfy a Dyson 
equation 

G a = + G^So-Go- = G CT + Go-So-G^.. (23) 

Summation over planar diagrams— yields the following 
expression for the self-energy matrix 



(24) 



Tr(Gi 2 ) 
Tr(G 21 ) 



® 1m, 



(25) 



in which G 12 and G 21 are off-diagonal blocks of the full 
Greens function 



Gn 



u a u 

r<2i r<22 



(26) 



The function £(z) is defined by the planar series, £(2) = 
Er=i w n z " _1 7 m which w n are symmetry coefficients 
given by 



(_l)n+i(2 n -2)! 



n\(n- ly.M 2 ™- 1 ' 
Evaluating the sum we obtain the closed expression 



(27) 



z 2 ((z 2 ) = i (V4z 2 + Af 2 - M 



Using l|28[) we may write Eq. 124|l in the form 

t a P a ± a +M± a = P a . 
To proceed, we introduce the variables 



ha = ^Tr(G* 2 ) 



0-6*1 



and 



N 



Tr(G 21 )=a0 2 



(28) 



(29) 



(30) 



(31) 



which implies the following simple form for the matrix 

Pa 



Pa = Na 







® 1m- 



(32) 



Using the ansatz 



(33) 



in Eq. I|29|l . we obtain the following coupled system 



(3 2 2 + 2p! = X 
(3 2 6 1 + 2f3 2 = 8 2 . 



(34) 



At this stage we need to further specify the nature of 
the contacts between the cavity and the leads. We as- 
sume the presence of asymmetric barriers with N equiva- 
lent transmitting channels. The adequate choices for the 
matrices T, L and R are therefore given by 



T = 



iy/Tl 







l;v = L 



(35) 
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and 



A. Chaotic Cavity with Symmetric Barriers 



R = 



vT=T7 

o yr 



(36) 



in which Ti and r 2 represent the transmission coeffi- 
cients of each channel in barriers 1 and 2 respectively. 
Combining these expressions with the ansatz l|33[) for the 
self-energy matrix, we may calculate explicitly the full 
Green's function, from which (after using Eq. wc 
obtain 



and 



aT 2 + 


(1 - r a )/? x 




1 - a0 2 T 2 
(1 


-/Wi- 


r 2 ) 


1 — afiiTi 




ro 


oTi + 


(i-ri)/3 2 




1 — afliTi 
(1 


-P1P2O-- 
- r 2 )/3 2 





+ 



(37) 



+ 



a/3 2 r 2 - /3i/3 2 (l - r 2 ) 



(38) 



Combining l|37|) and (|38|l with 134|) , we obtain the follow- 
ing non-linear algebraic system 

a(i-r 1 )ry? 1 /3 2 3 + 
[(a 2 r x r 2 + {2v l - 1) r 2 - r x ) X - a (i + r x ) r 2 ] /3|+ 
[(ri - 2r!r 2 ) aft + a 2 rir 2 + r : + r 2 ] 2 - a v l = o, 

(39) 

and 

a (i-r 2 )r 1 ^+ 
[(a 2 r!r 2 + (2r 2 - 1) r x - r 2 ) & - a (l + r 2 ) /3 2 + 

[(r 2 - 21^2) a(3 2 + a 2 rir 2 + 1\ + r 2 ] ft - a r 2 = o. 



(40) 



Finally, inserting the ansatz, Eq. (|24H . into (|18|) yields 

oT 2 /3 2 



and 



h{z) = a z N 



f 2 (z)=a 2 N 



1 - 



1 



1 - (1 - r 2 ) ftft 



1 - (1 - r x ) fafa 



(41) 



(42) 



Equations l|39"j) , (gUJ) , l|4*2|) and (|4*T|) are the central results 
of this section. Together with Eq (|16fl they represent 
the complete solution of the problem. It extends the 
calculations of Brouwer and BeenakkerJ 1 - to include the 
case of asymmetric barriers (although in their solution 
of the symmetric case the channels in the barriers were 
considered non-equivalent). Before presenting a complete 
analysis of the generic situation, we shall discuss below 
some important particular cases. 



The case of symmetric barriers is described by setting 
T x = T = T 2 and [3 l = = 2 in Eqs. GJj), £P, (£E1> 
and l|41|) . so that we find 



(a0 2 -2(3 + a)((l - T)(3 2 + aT(3 - 1) = 0, 



(43) 



and 



ji(z) = a 2 7V 



1 - 



1 - (1 - r) p 2 



= f 2 (z). (44) 



The physical root of Eq. (f^ is given by (3 = (1 — 
Vl — a 2 ) /a = \fz— v 1 z — 1, which yields fi(z) — f(z) — 
/ 2 (z), where 



f{z) _ 2 (y/i - y^T) (1 - r) + r/7J^T 
iv iz{yfl- vF=~i)(i-r) + iyi 

Inserting l|45|l into (|16[1 we find the average density 

iV r(2-r) 



P(r) = 



7T (r 2 -4rr + 4r) v /r(l - r) ' 



(45) 



(46) 



in agreement with Ref. Il3j The solution obtained by 
Brouwer and Beenakke»ii is, in fact, a generalization of 
Eq. Q46JI for non-equivalent channels. It reads 



JV 



r n (2 - r T . 



(47) 



^ 7T(T 2 n - 4r„r + 4r) v /r(l-r) 
which, of course, reproduces l|46[) when r ra = T for all n. 

B. Chaotic Cavity with Two Tunnel Junctions 

This system is described by applying the conditions 
Ti,r 2 < 1 in Eqs. HP, gUJ), g2J and 62). We get 



/3lar 2 - ( s 2 (ri + r 2 ) + ari = o 
/3 2 ar 1 -^i(r 1 + r 2 ) + ar 2 = o, 



together with 

h{z) =a 2 N 

and 

f 2 (z) = a 2 N 



aT 1 1 



-1 -1 



1 - 0x02 



aT 2 (3 2 



1 ~ 0102 



(48) 



(49) 



(50) 



From the physical roots of (|48|l . we obtain fi(z) = f(z) 
/ 2 (z), where 



N 

z 



rir 2 



(Ti + r 2 )^/z(z-T ) 



(51) 
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in which r = 4T 1 T 2 /(T 1 + T 2 ) 2 . Inserting (J5TJ into lO 
we obtain 



ATxr 2 



1 



7r(r 1 + r 2 ) r 3/2 v /^~7)' 



(52) 



in agreement with Refs. [l2 and Ha. 

Using 1321) we can compute the average value of several 
physical observables. Of particular interest are the first 
three cumulants of the FCS, defined in Eqs. 10, © and 
©. We find 



{9)i 



N- 



11 2 



for the average conductance, 



r 2 + r 2 



MJ T j (Ti + r 3 )2' 

for the Fano factor, and 

( K )\ rf - 2r?r 2 + 6r 2 r 2 - 2^ + r| 



(p) 



TJ 



(r 1 + r 2 ) 2 (r 2 + r 2 ) 



(53) 



(54) 



(55) 



for the ratio between the average third cumulant and the 
average shot-noise power. The subscript TJ stands for 
Tunnel Junction . Equations i|53|) , l|54|) and (|55|l are well 
known in the literature^. 



III. GENERAL ANALYSIS AND COMPARISON 
WITH CIRCUIT THEORY 

In this section we compare the predictions of the above 
diagrammatic approach with those of the extended ver- 
sion of circuit theory. The exact agreement that we found 
for various quantities strongly suggests the full equiva- 
lence of these semiclassical techniques in the entire range 
of input parameters. 



A. Circuit Theory 

Circuit theory was invented by NazaroAii^ and repre- 
sents a very powerful tool for computing ensemble av- 
erages of quantum chaotic systems in the semiclassical 
regime. It consists of a finite element approach in which 
the spatial support of the system is partitioned into a 
network, containing edges (or connectors) and vertices 
(nodes or terminals). In its simplest version, there are 
only two terminals and the system is subject to a ficti- 
tious pseudo-potential, with fixed values at the termi- 
nals and unknown values at the internal nodes. The basic 
principles are: 1) a general law of pseudo-current- voltage 
(I — V) characteristics for a connector (i, j) subject to a 
pseudo-potential drop 



J«(A*«) = / dr 



1 -Tsin 2 (A$y/2)' 



(56) 



where 2y(A$y) is a pseudo-current and pij(r) — 
(TV(5(r — 1S the average transmission eigenvalue 

density of the connector. We defined tij as the sample 
specific transmission matrix of the connector. 2) a gen- 
eralized Kirchhoff law for pseudo-current conservation at 
each node. 

The power of this approach depends crucially on the 
appropriate choice of the connectors, which in turn 
depends on our ability to calculate the correspondent 
pseudo I — V characteristics. This is the point where the 
extended approach, put forward in Ref. 0^, differs from 
Nazarov's original formulation. In Ref. |lj the function 
Iij(A$>ij) was calculated directly from the supersymmet- 
ric non-linear cr-model. The advantage of this procedure 
is the possibility to take full account of quantum effects in 
the description, whenever necessary, such as near barriers 
and interfaces, without the need to introduce additional 
assumptions for performing the ensemble averages. 

As an example, consider the system studied in previ- 
ous section, i.e. a double barrier chaotic billiard with 
N equivalent channels at each contact. In the extended 
circuit theory, we end up with only two equations 13 



i(0) = i 1 (4>-o) = i2(o), 



where 



W) = 



1 



NTj 8w(<f>) 
-Tj sin 2 (0/2) 



(57) 



(58) 



is the pseudo I — V characteristics of barrier j, inter- 
preted as a connector. We remark that these equations 
contain all quantum information that remain relevant in 
the semiclassical regime after ensemble averaging. They 
are therefore a direct consequence of the maximum en- 
tropy principle. In Ref. it was found to be convenient 
to introduce the following modified pseudo-current 



K(x) = -I(-2ix), 

which yields the conservation law 

K(x)=K 1 (x-y)=K 2 {y), 

where 



1 



1 



tanh(x + — otj) + tanh(a; — —a 



(59) 



(60) 



(61) 



in which we introduced the constants, ctj, via the relation 
Tj = sech 2 (a J 72). Inserting (JSTJ into jjjOJ yields 

tanh(a; — y + ~^ a i) + tanh(x — y — 2 a i) = 



tanh(y + ^a 2 ) + tanh(y - ^a 2 ) 

which after using the trigonometric identity 

tanh x — tanh y 



tanh(x — y) 



1 — tanh x tanh y ' 



(62) 



(63) 
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yields the following polynomial equation for the variable 
£ = tanh y 

[Ti(l - r 2 )tanhx]£ 4 - (3rir 2 tanhx)£ 2 + 

[(Tir 2 + r 2 - r x ) tanh 2 x + 2r!r 2 - r x - r 2 ]£ 3 + 
[(r x r 2 + r x - r a ) tanh 2 x + r x + r 2 ]£ 

-ritanha; = 0. (64) 
This equation must be supplemented by 

m = ^^Ue - (65) 

Equations <(64(1 and 1(65(1 were first presented in Ref. 
and are the circuit theory equations for the double barrier 
chaotic billiard. In the next subsection we shall compare 
them with Eqs. 1(23), 601, (32) an d J3H, obtained from 
the diagrammatic method. 



B. Comparison between Circuit Theory and the 
Diagrammatic Approach 

In order to facilitate comparison, let us first estab- 
lish a direct relation between the generating function 
fj(z) = TrFj(z) of the diagrammatic technique and the 
modified pseudo-current K(x) of circuit theory. From 
((13(1 we obtain 

oo 

F 1 (z) = Y,(c 1 (tfC 2 SC 1 ) n ) — T (66) 



n=0 



and 



Inserting J2J and (|15|l into the above equations we get 



oc 1 



n=0 



Z n+1 \ (tU) n 



(68) 



and 



1 /(ttt)" 



n=0 



^ 



Performing the trace, we obtain fi(z) — f(z) — f 2 (z), 
where 

z 2 f(z) =Nz + h(l/z), (69) 
and we introduced the function h(z), defined as 

From (|56|l . one realizes that 

7(0) = /i(sin 2 (< ? ! ) /2))sin0, (71) 



which when combined with (|69|l yields 



f(z) = -(N + 



2y/7^T 



sin(<A/2)=l/Vi 



(72) 



Using the relation between K(x) and /(</>), Eq. 1(59(1 . we 
may rewrite the above equation as 



/(z) = -[ N-K(x) tanh a; 



(73) 



sinh x—\j s/—z 



Let us consider two simple applications of Eq. 1(73(1 . For 
a chaotic cavity with symmetric barriers, Eqs. 164|) and 
(E3 yield 



isT(a5) = 



AT sinh x 



2-r + Tcosha;' 
Inserting JJ3J mto G3 w e find 



(2-r)y^T + ryi 
I V^T" [(2 - r) + rVT^T] 



(74) 



(75) 



which can be easily shown to agree with l|45|l . For a 
chaotic cavity with two tunnel junctions, we obtain from 
Eqs. 1)64(1 and 1(65(1 the following expression 



K(x) = 



AT^ sinh2x 



2y/Tl + T| + 2rir 2 cosh 2x 
Inserting 1)76(1 into ((73(1 we obtain 



(76) 



/(*) = 



N 



1^ 2 



(ri + r 2 )^(^-r ) 



(77) 



in agreement with Eq. 1(51(1 . 

We can also invert Eq. ((73(1 to obtain K(x) from the 
diagrammatic approach 



K(x) = coth(x) TV - z/(z) 



(78) 



z — — sm 



h- 2 (x) 



In particular, we can calculate the average of the cumu- 
lants of the FCS using the formula 



h k+1 = (Tr{(ttf) k+1 }) 



(-l) fe 2 fc d k H{x) 



ft! d(cosh(2x)) fe 



x=0 



(79) 

where H(x) = 2K (x)/ sinh(2x). The first three cumu- 
lants are then given by 

{g) = h u (p) = h 1 -h 2 , (k) =h 1 -3h 2 + 2h 3 . (80) 

We remark that the function if (x) can also be obtained 
directly from f(z) via the relation 



H(x) =z(zf(z)-N 



(81) 



z—— sinh 2 (a:) 
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Motivated by strong numerical evidences^*, our basic 
conjecture is that one obtains the same function H(x) 
either from Eqs. l)31fy. jiffy. l(l2fy and (gTJ or from Eqs. 
(|64l) and (|65|l . Although explicit analytic evaluation of 
-ff(x) in both approaches, for general Ti and r 2 , is too 
cumbersome (albeit possible in principle), we can still 
make some analytical progress in verifying our conjec- 
ture by expanding H(x) in powers of x, and evaluating 
the coefficients as explicit functions of Ti and r 2 . This 
procedure yields, in both approaches, the following closed 
expressions 



(9) = N 



rir 2 
ri + i 



= (9) 



TJ 



for the conductance, 



M 
(9) 



rir 2 



(9)J T j 



for the Fano factor, and 



(P) 



2r x r 2 



(p) 



(82) 



(83) 



(84) 



TJ 



for the ratio between the average third cumulant, (k), 
and the average shot-noise power, (p). We also found 
agreement for the fourth cumulant. Note that for tunnel 
junctions, we have Ti,r 2 <C 1, and 1(52 )1 .1(53 )1 and ((51)) 
reduce to IpHfy. ((51) and J55J), as expected. The above 
expressions together with the forth cumulant are in com- 
plete agreement with the semiclassical cascade approach 
presented in Ref. They represent strong analytical 
evidences for the mathematical equivalence between Eqs. 
((23), 601, and (gTJ of the diagrammatic technique 
and Eqs. H64|) and l|65|) of circuit theory, in the descrip- 
tion of the asymptotic semiclassical domain of the Pois- 
son kernel distribution. 



ysis of the unitary group and circuit theory, in the de- 
scription of quantum transport through double-barrier 
chaotic billiards. The problem was reduced to a com- 
parison between a pair of coupled non-linear algebraic 
equations (diagrammatic technique) and a polynomial 
equation of fourth order (circuit theory). Exact agree- 
ment was found for a variety of quantities, such as the 
first four average cumulants of the full counting statis- 
tics and the average transmission eigenvalue density for 
symmetric barriers and tunnel junctions. The complete 
equivalence of these approaches is a non-trivial result, 
because the semiclassical concatenation principle, used 
to derive circuit theory equations, has no obvious coun- 
terpart in the diagrammatic method and leads to a sub- 
stantial algebraic simplification of the whole problem. It 
is quite amusing to observe that circuit theory turns out 
to yield a irreducible description that intuitively satisfies 
the "Occam's razor" principle of descriptive simplicity 19 . 

An interesting consequence of our result would be the 
extension of circuit theory to deal with quantum correc- 
tions, which can be systematically treated in both the 
supersymmctric non-linear cr-model and in the diagram- 
matic technique. From a broader perspective, we expect 
our result to help establishing a direct connection be- 
tween several recent independent developments of both 
circuit theory and the diagrammatic S'-matrix approach, 
such as those in magnetoelectronics^ with obvious tech- 
nical advantages. 



IV. SUMMARY AND PERSPECTIVES 

We presented a detailed comparison between two well 
known semiclassical approaches, the diagrammatic anal- 
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